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Chemical adsorption of atomic hydrogen on a negatively charged single layer graphene sheet has been ana- 
lyzed with ab-initio Density Functional Theory calculations. We have simulated both finite clusters and infinite 
periodic systems to investigate the effect of different ingredients of the theory, e.g. exchange and correlation po- 
tentials, basis sets, etc. Hydrogen's electron affinity dominates the energetic balance in the charged systems and 
the extra electron is predominantly attracted to a region nearby the chemisorbed atom. The main consequences 
are: (i) the cancellation of the unpaired spin resulting in a singlet ground-state, and (ii) a stronger interaction 
between hydrogen and the graphene sheet. 
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I. INTRODUCTION 

Planar two-dimensional graphene has been considered to 
be a very promising new material since its preparation by 
Novoselov et alJ^ and Berger et al^ in 2004. Recent interest 
is focused on the appearance of magnetism around point de- 
fects in graphene 3 and the possibility of hydrogen storag 
A global understanding of the origin of magnetism in finite 
graphene systems is provided by Lieb's theorem on bipartite 
lattices 8 . Any unbalance between the numbers of sites be- 
longing to each of the two sublattices gives rise to a mag- 
netic groundstate 9 . This result rests on the validity of a sim- 
ple Hubbard Hamiltonian which certainly works for the semi- 
quantitative description of the spin states of some polycyclic 
aromatic hydrocarbon molecules but not for the correspond- 
ing charged states 10 . The saturation of a carbon 7r-electron by 
hydrogen is one of the simplest ways to change the balance 
between sublattice sites and, consequently, it produces spin 
polarization in the neighborhood. This provides an interest- 
ing link between hydrogenation and spin production that mo- 
tivates this work. We shall demonstrate that the 1/2 spin orig- 
inated by the presence of an isolated hydrogen atom on top of 
a carbon atom belonging to planar graphene can be quenched 
if an extra electron bounds to the defect. Numerical results 
obtained by ab initio open-shell Density Functional Methods 
(DFT) suggest that the complex defect is energetically favor- 
able. A similar spin quenching phenomenon was discussed 
some years ago by Duplock, Scheffler, and Lindan for hydro- 
gen near a Stone-Wales defect 11 . These authors have shown 
that the spin polarized groundstate around chemisorbed hy- 
drogen disappears in the presence of a Stone- Wales defect. In 
their interpretation, this result is probably due to the strong 
destruction of alternation near the defect that eliminates the 
tendency to antiferromagnetic order. In our case, however, we 
argue that the mere flow of charge is enough to heal the lo- 
cal unbalance between sublattices due to the existence of one 
chemisorbed hydrogen on an otherwise ideal graphene mate- 
rial. 

The rest of the paper is organized as follows: The ab initio 
methods used in this work are presented in section II followed 
by the discussion of our main numerical results. Final section 
III just remarks our main message. 



II. AB INITIO CALCULATIONS: METHODS AND 
RESULTS 

Two planar carbon clusters with the structure of graphene 
have been chosen to calculate the energetics of hydrogen ab- 
sorption both at the neutral state and when the system is elec- 
tronically charged with an extra electron (anions). The first 
system is the well known coronene polycyclic aromatic hy- 
drocarbon (PAH) represented in Fig.[T]while the second one is 
a larger PAH obtained from coronene adding an extra shell of 
benzene rings. This system is referred as supercoronene in the 
literature. We notice that this larger PAH molecule has not yet 
been synthesized but nonetheless it provides a good theoreti- 
cal benchmark for our purposes (Fig. [2] schematically shows 
one hydrogen chemisorbed on supercoronene). The ground- 
state for both PAH's does not show spin polarization (total 
spin is zero). This is an important difference with our previ- 
ous study of hydrogen chemisorption on graphene where spin 
one-half clusters were used to preserve the point symmetry 
and facilitate the computational effort 5 . In the present work, 
however, we focus on the description of spin polarization and 
we must start with an unpolarized cluster to accurately simu- 
late the graphene layer. Point symmetry is lost and computa- 
tional load is larger. Nevertheless, we have checked that both 
structural and energetic results for chemisorbed H coincide 
with the values given in our previous worP. 

Quantum chemistry calculations have been done using the 
GAMESS prograrrPl Several sets of gaussian basis functions 
have been employed. Depending on the computational ef- 
fort we have been able to assess the convergence of numer- 
ical results for some cases by comparing results obtained with 
bases of different sizes. For the largest systems, however, 
we have been forced to choose a minimal basis and conver- 
gence could only be assessed by reference to the smaller clus- 
ters. Specifically, we have always started by trying the so- 
called MIDI basisP, and when possible we have moved to 
a correlation consistent basis referred in the literature as cc- 
pVTzP (CCT within GAMESS and this paper) and a DFT 
adapted hierarchy of basis called PCn where n indicates the 
level of polarization^. Our best results correspond to the 
larger PCn basis that we have been able to use in each case. 
Unrestricted Hartree Fock (UHF) calculations have been in- 
tentionally avoided because total spin of the wavefunction is 
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TABLE I: Total energies (in hartree) for atomic hydrogen and its anion, for coronene (C24H12) and coronene anion (C24Hi2)~ both ideal 
and with a chemisorbed H atom and for a larger cluster (C54Hig) sometimes called supercoronene and its anion both planar and deformed by 
the presence of a chemisorbed H atom. Since results are given for two or more gaussian basis sets some rough estimation of error is possible. 



Cluster 


Total Spin 


Energy(MIDI) 


Energy(CCT) 


Energy(PC2) 


Energy(PC3) 
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1 
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-0.4953 


-0.4988 


-0.4990 


-0.4991 
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-2056.1072 




-2068.9095 
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-2056.1868 




-2068.9919 





undefined in those cases. Therefore, HF calculations for an 
odd total number of electrons have been performed using the 
Restricted-Open-Shell variant. All results presented in this 
work for clusters have been obtained using the Becke-Lee- 
Yang-Parr hybrid density functional 

RB3LYlfl6] 

Our choice of finite clusters of various sizes poses the ques- 
tion of to what extent results are affected by the particular 
boundary conditions imposed to solve the quantum problem 
for electrons. Therefore, we check by comparing with calcu- 
lations performed on extended periodic models using periodic 
boundary conditions and a plane-waves basis. This model is 
set up so a single H atom is adsorbed on a 4 x 4 graphene su- 
percell including 32 carbon atoms on a honeycomb lattice (see 
Fig.|6] a = b = 9.84 A, c = 23.4 k,a = fi = 90°, 7 = 60°). 
We use ultrasoft pseudopotentials^, an energy cutoff of 310 
eV, and a Monkhrost-Pack mesh of 3 x 3 x Actual calcu- 
lations are performed with the CASTEP program allowing for 
spin polarization of the different electronic bands 19 20 . To de- 
scribe the exchange and correlation potential we use the local 
density approximation 21 . Forces and stresses on the system 
are converged to the usual thresholds (9 x 10~ 3 eV/A and 0.01 
GPa) and the total energy is minimized for the different sys- 
tems. This procedure cannot provide an accurate description 
for the electron affinity energy since the extra electron in the 
supercell is intentionally neutralized with a uniform positive 
background to subtract infinite contributions in the periodic 
system. This uniform background, however, has little effect 
on the spatial distribution of the electronic and spin densities, 
that can be analyzed with reasonable confidence. 

Table I compiles the bulk of our quantum-chemistry re- 
sults. Notice that together with the total energies needed to get 
chemisorption energies for hydrogen on graphene, we have 
computed the energies corresponding to systems charged with 
an extra electron (anions in the molecular case). Total spin 
of the groundstate is given by the second column of the ta- 
ble. It can be seen that the S=0 value of the neutral molecules 
is recovered by the anions of the hydrogenated cases. Table 
II gives the electron affinities that are obtained from the re- 
sults shown in Table l 22 . Nice convergence to the experimen- 
tal electron affinity of hydrogen is observed in the first entry of 
the table. Results for larger clusters are limited by our compu- 



TABLE II: Electron affinities of hydrogen, coronene, monohy- 
drogenated coronene, supercoronene and monohydrogenated super- 
coronene obtained from the results compiled in Table I. Energies are 
now given in eV. 

Cluster Energy(CCT) Energy(PC2) Energy(PC3) Experimental 



H 


0.13 


0.51 


0.72 


0.75" 


C24H12 


0.32 


0.39 


0.40 


0.47* 


C24H13 


1.48 


1.57 






C54H18 


1.07 r 


1.35 






C54H19 


2. 17'' 


2.24 







"See, for example, http://www.chemicool.com/elements/hydrogen.html 
*The electron affinities of coronene and coronene dimer have been measured 

by Duncan et at as reported in Ref. 423t . 
c This result corresponds to the small MIDI basis type of calculation that is 

described in Ref.JTJJ. 



tational means but nonetheless our results at the PC2 level are 
good enough to support our conclusions. A word of caution 
is in order here: although these results for the spin seem to 
fit nicely within a simple electron-counting scheme (i.e., zero 
spin for even number of electrons and net spin for odd number 
of electrons) this is not always true. In particular, we recall the 
case where two hydrogen atoms are adsorbed on the graphene 
layer: while adsorption of the two hydrogens in next-neighbor 
positions results in a ground state with no net spin, adsorption 
in next to next-neighbor sites results in a ground state with 
a net S — 2 fj,B- This is related to the fact that graphene is 
a bipartite lattice, and it is in accordance to Lieb's theorem 8 , 
showing the limitations of simple electron-counting rules. 

Let us briefly discuss the results given in Table II. The elec- 
tron affinity of C24H43 cluster, i.e., the cluster with one hydro- 
gen chemisorbed on top of one of the six carbon atoms on the 
inner ring of coronene (see, Fig. 1) is 1.18 eV larger than the 
electron affinity of coronene. Also, the value for C54H19, that 
is, H on top of supercoronene is 0.89 eV larger than the super- 
coronene electron affinity. Although, only two cluster sizes 
have been studied, it seems that the difference is approaching 
a limiting value close to the electron affinity of free hydro- 
gen, that is, close to 0.75 eV. If this were the case, the extra 
charge would be attracted by hydrogen with a similar strength 
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as in free space and the screening by the rest of 7r-electrons 
of graphene would remain unnoticed. Nevertheless, we show 
later that only part of the extra charge remains close to the de- 
fect. Therefore, we assume that for an infinite system only a 
fraction of 0.75 eV proportional to the localized charge would 
remain. 

There is an alternative elaboration of the results given in 
Table I focusing on the variation of the binding of a hydro- 
gen atom on top of a charged surface compared to the binding 
by the neutral one. From this point of view, H binding en- 
ergy increases from 0.59 eV to 1.77 eV on coronene and from 
0.36 eV to 1.25 eV in supercoronene using PC2 values in Ta- 
ble H This means that the charged systems bind hydrogen 
about one eV stronger than the neutral ones. From this per- 
spective, the larger values of the electron affinity obtained for 
hydrogenated clusters can be asigned to a stronger hydrogen 
binding to graphene. 

The relative facility to move charge across the overall sys- 
tem implied by the semi-metallic character of graphene and 
our results in Table II point towards the formation of a com- 
plex point defect with an extra electron in the neighborhood 
of the chemisorbed hydrogen atom. This picture is further 
supported by the spatial distribution of this extra electron in 
the studied clusters as it is shown in Figs. [3] and [4] Charge 
densities obtained with PC2 gaussian basis for coronene an- 
ion and neutral coronene have been subtracted in Fig. fus- 
ing the MOLD EN package 25 . The same difference for super- 
coronene is given in Fig. |4] In both cases, the charge around 
"on top" H is similar to the extra charge of the H anion. On 
the other hand, the spreading positive and negative densities 
are similar but not equal in coronene and supercoronene. A 
closer inspection reveals that the extra electron is occupying 
the partially occupied HOMO level of the corresponding neu- 
tral clusters. Fig. [5] gives a picture of the HOMO orbital of 
H on supercoronene that nicely explains the charge difference 
previously shown in Fig. |4] It is interesting to notice, how- 
ever, that although the extra charge clusters around the adsor- 
bate, it is partly delocalized as it is expected from quantum 
mechanics principles. In fact, from independent tight-binding 
calculations in periodic systems we don't find a true expo- 
nential localization around the defect (see Appendix). On the 
contrary, a percentage of the charge is extended all over the 
system (e.g., see Fig. [4]). However, as our detailed quantum- 
mechanical calculations show, the net effect of the localized 
part is enough to quench the spin. 

Our ab initio results on periodic extended systems fully 
support the interpretation given in the previous paragraph. We 
observe in Fig.[6]how the extra charge is accumulated around 
the adsorbed H. In this case, bonding charges for the neutral 
and charged supercells are depicted in the upper and lower 
panels of the figure, respectively. Bonding charges are de- 
fined as density charge differences between the whole system 
and conveniently defined fragments. In our case, hydrogen 
atom is one fragment while the graphene 4x4 supercell is the 
second one. The upper difference integrates to one electronic 
charge since the whole system is charged while the fragments 
are neutral. We notice that a small part of the extra electron is 
on carbon atoms while an important part attaches to hydrogen 



(the violet negative density in the lower panel does not appears 
in the upper panel meaning a positive contribution to the extra 
charge in the upper panel). The same picture is extracted from 
a Mulliken analysis of populations around different atoms. In 
the neutral system, charge flows upon adsorption from hydro- 
gen to graphene, so approximately —0.63 e is located around 
hydrogen while the transferred charge resides mostly around 
the closest carbon behind hydrogen (—0.33 e). On the other 
hand, in the charged system we find — 1.42 e around hydrogen, 
while the carbon behind keeps nearly the same occupation 
(—0.34 e) and the rest of the charge is distributed over nearest 
neighbors and next-nearest neighbors. Therefore, about 80% 
of the extra electron is located near the chemisorbed hydro- 
gen. Along the same line, the H-C bond population analysis 
is about three times larger for the anion, although the bond 
length is nearly not affected. Finally, integration of the spin 
density and the absolute value of the spin density over the 
simulation cell give further support for this picture. In the 
neutral system these values amount to 0.4 and 0.5 \lb respec- 
tively, while in the charged one decrease to values of amount 
to 4 x 10~ 6 and 5 x 10~ 4 /is respectively. Therefore, the ac- 
cumulation of charge around the adsorbed H and the C nearest 
and next nearest neighbors plays the role to cancel the extra 
spin polarization brought by the adsorption of H on the clean 
graphene layer in accordance with the results obtained on fi- 
nite clusters. 

There is a subtle chemical argument that helps the under- 
standing of our numerical results. In Ref. (l26t . trivalent car- 
bon atoms with an unpaired electron were unraveled in the 
studied geometry of carbon tetrapod. This carbon radicals 
were stabilized by steric protection giving rise to unpaired lo- 
calized electrons that polarize the carbon neighborhood and 
explain the appearance of magnetism in pure organic sys- 
tems. The original paradigm is tryphenilmethyl, synthesized 
by Gomberg in 190C 27 , where a trivalent carbon atom is sta- 
bilized by three bonds to benzene rings impeding the reaction 
with a similar molecule. Nevertheless, the anion of tryphenil- 
methyl reacts with a proton to form a strong C-H bond (heat 
of reaction amounts to 15.65 eV per molecule): 

C19H15 + H + — > CigHi 6 , 

producing a neutral non-magnetic molecule resembling the 
clusters that we have studied here^. Both the number of H 
atoms and electrons are even allowing an easy shell closing 
and stabilization of the resulting molecule. We can adapt the 
underlying chemistry of these phenomena to the binding of H 
on graphene using the following argument: it can be thought 
that when H forms a covalent bond with a C atom of graphene 
the two binding electrons are paired but additionally and due 
to the particular topology of graphene lattice one 7r-electron 
becomes unpaired and, therefore, spin polarized. Since there 
is no steric protection around the defect (a kind of radical) 
any more or less free electron of the system will flow into the 
defect to restore the equilibrium between sublattices. Conse- 
quently, spin polarization around chemisorbed hydrogen dis- 
appears. This qualitative argument is fully supported by our 
total energy calculations showing a gain in potential energy 
following the spin neutralization (remember that the binding 
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energy of H increases about one eV for the charged system). 



m. CONCLUDING REMARKS 

From a detailed analysis based in first-principles DFT cal- 
culations we find that in the presence of an extra electron 
chemisorbed H plays to keep most of the extra charge in its 
vicinity. The electron affinities computed on finite cluster 
models seem to converge to a value that is somewhat smaller 
than the free atomic hydrogen value of 0.75 eV. Our calcula- 
tions suggest that being graphene a semimetal with zero den- 
sity of states at the Fermi energy, screening of Coulomb inter- 
actions by the 7r-electrons liquid is weak and allows the elec- 
tron flow to sites where H is chemisorbed, forming a complex 
point defect. Accumulation of extra charge around the defect, 
otherwise giving rise to spin polarization, works to quench it. 
We suggest that this physical effect is behind the difficulties 
to observe magnetism in graphene derived systems. 
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Appendix: A tight-binding calculation 

A tight-binding model including carbon 7r-orbitals and hy- 
drogen Is orbital allows a straightforward although approxi- 
mate solution of the topic covered in this paper. Since larger 
sizes can be studied, qualitative conclusions on the scaling be- 
havior of the charged system can be obtained. This calculation 
fully supports our arguments based on the computationally de- 
manding ab initio study done in the main part of the paper for 
smaller systems. Although total energies are not attainable by 
the method, both the density of states and the charge distri- 
bution of an extra electron suggest the stability of a charged 
hydrogen atom bonded to graphene. 

The p^-p z hopping term of the Hamiltonian is -2.71 eV. 
Comparing the ionization energy H (13.60 eV) to the value 
for the methyl radical (9.84 eV) we take the H level relative to 
the 7r one as -3.76 eV. This is justified because of the similarity 
between graphene and methyl radical p z orbital. For the C-H 
coupling we have taken -5.39 eV based on a simple scaling 
argument. A standard supercell procedure and a carefully BZ 
integration have been used to get accurate density of states 
and local charge magnitudes. Although very large supercells 
have been explored (50 x 50), we will give here results for a 
10 x 10 system because it shows almost converged results that 
are better depicted. 

Fig.([7]i shows the total density of states per spin of the su- 
percell. It integrates to 402 (total number of states of the 
cell) and shows two salient features: (i) a deep level below 
the valence band describing the bonding C-H orbital, and (ii) 
strongly perturbed values below Dirac's Point energy (0 in the 



clean system). Nevertheless, the most interesting feature rela- 
tive to this density of states (DOS) is that it integrates to 202 
up to the Dirac energy (E=0), i.e., the complete occupation 
of the valence band describes one extra electron above half- 
filling which is only 201 (200 electrons corresponding to the 
200 7r-orbitals of C atoms plus one electron contributed by 
H). This means that when the Fermi level coincides with the 
Dirac point of graphene, a region with a quimisorbed H atom 
becomes charged by an extra electron. Since the full occu- 
pancy of the valence band can be thought as closing electronic 
shells, some stabilization of the system can be inferred from 
this property of the DOS. 

The ideas suggested by the DOS are confirmed by the spa- 
tial distribution of the electronic charge induced by the addi- 
tion of an extra electron to the system. Fig. ([8]) (upper panel) 
has been obtained subtracting the background charge of one 
electron per site from a 10 x 10 supercell with an additional 
electron. The resulting distribution can be compared with the 
corresponding variations of Mulliken populations of our ab 
initio results (lower panel). It shows that an important amount 
of charge directly resides on the hydrogen atom although the 
underlying carbon atom somewhat reduces the net value on 
the binding site in the empirical model results. Yet another 
important part of the electron lies close to this site whereas 
the rest is spread over the whole supercell. The most intrigu- 
ing characteristic of the charge distribution is its dual nature 
of localized around the defect and extended over the whole 
system 29 . In any case, a tight binding model with a minimum 
number of orbitals reinforces our idea signaling the tendency 
of electrons to remain in the neighborhood of bounded hydro- 
gen. 
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FIG. 1: (color online) Upper panel: schematic representation of a 
coronene molecule showing a honeycomb lattice inner structure sat- 
urated in the boundary by hydrogen atoms so the coordination of 
carbon atoms is preserved over the whole system. Lower panel: An 
extra hydrogen atom is chemisorbed on top of a carbon atom belong- 
ing to the inner ring of coronene. Although the sp 2 to sp 3 reconstruc- 
tion is only faintly visible in the figure, the C-H bonding distance and 
the details of the upward relaxation of C neighboring atoms coincide 
with those given in Ref.([5j. 



FIG. 2: (color online) Hydrogen chemisorbed on supercoronene 

(C54H18) 




FIG. 3: (color online) Charge difference between the anion of 
coronene and the neutral molecule. Density iso-contours of ±0.002 
e/A 3 are represented. 



FIG. 4: (color online) Charge difference between the anion of super- 
coronene and the neutral molecule. Density isocontours of ±0.002 
e/A 3 are represented. 




FIG. 5: (color online) HOMO orbital of supercoronene. Isocontours 
of ±0.05 A are represented. 



FIG. 6: (color online) A 4 x 4 graphene supercell with a single ad- 
sorbed H atom showing the bonding charge density first for the sys- 
tem with one extra electron (upper panel) and, second, for the neutral 
system (lower panel). Density isocontours of ±0.04 e/A 3 are given 
(blue and violet respectively). The accumulation of charge near H is 
seen to be the origin for the cancelation of the extra spin. 
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FIG. 7: Total density of states of a 10 x 10 supercell calculation of a 
graphene layer containing one hydrogen atom forming a strong cova- 
lent bond with one of the carbon atoms of the layer. Its integration up 
the the Dirac point at E=0 gives one extra electron above half-filling. 



FIG. 8: (color online) Spatial distribution of the extra electron charg- 
ing the hydrogenated supercell. The sphere volume is proportional to 
the extra number of electrons of the corresponding atom, red if pos- 
itive and blue if negative. The relevant part of the tight-binding re- 
sults obtained for a 10 x 10 supercell calculation is given in the upper 
panel whereas ab initio results for the 4x4 supercell are shown in the 
lower panel. Last values are based on Mulliken populations. Overall 
trends coincide but details differ: both bonded H and C atoms get 
extra population when analyzed by the Materials Studio (CASTEP) 
package but not for the semiempirical model, for example. Although 
this representation of the extra charge cannot be directly compared 
with isocountours shown in Figs. ([3} and {4}, the overall distribution 
looks similar. 



